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f *^ 1 We study in this paper the sine-Gordon model using functional Renormalization Group 

(fRG) at Local Potential Approximation (LPA) using different RG schemes. In d = 2, 
using Wegner-Houghton RG we demonstrate that the location of the phase boundary 
is entirely driven by the relative position to the Coleman fixed point even for strongly 
coupled bare theories. We show the existence of a set of IR fixed points in the broken 
phase that are reached independently of the bare coupling. The bad convergence of the 
Fourier series in the broken phase is discussed and we demonstrate that these fixed- 
points can be found only using a global resolution of the effective potential. We then 
1 introduce the methodology for the use of Average action method where the regulator 

$_l ' breaks periodicity and show that it provides the same conclusions for various regulators. 

The behavior of the model is then discussed in d ^ 2 and the absence of the previous 
fixed points is interpreted. 
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1. Introduction 

The sine-Gordon model and its phase diagram has received a considerable interest 
these last decades due to its unique properties. In the euclidean spacetime, its 
Lagrangian reads : 

C=^(d^) 2 +ucos(f3<f > ) (1) 

Most of our knowledge has been established in two dimensions. It is has been 
shown perturbatively -for u « 1- to exhibit two phases separated by a Kosterlitz- 
Thouless transition 33 . The spontaneously broken periodicity phase for f3 < fi c with 
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asymptotic freedom and a non-renormalizable one for (3 > f3 c separated by the 
Coleman frequency (3 C = \/&k 20 . 

The spontaneously broken phase is of major interest as it has been demonstrated 
to be equivalent to the neutral sector of the massive Thirring model at zero 20 > 41 and 
finite temperature 23 . At the partition function level, the sine-Gordon model can 
also be mapped onto the neutral Coulomb gas 61 ' 8 . The sine-Gordon model has also 
been demonstrated to be dual to the XY model in the Villain approximation 63 for 
vanishing vortex fugacity 29 . And finally its ^-matrix has been exactly computed 
in 9 as well as the ^-matrix for its solitons 68 . 

In higher dimensions, the mapping at the partition function level between sine- 
Gordon and the Coulomb gas still holds as well as its relation to the non-linear sigma 
model in the small explicit breaking limit 61 but little is known on these models, 
especially in the IR of their broken phase where non-perturbative techniques use is 
mandatory 32 . 

In addition to the previous relations with other models, the sine-Gordon model 
is also the simplest non-trivial theory with a field compact variable, a feature shared 
with non-abelian gauge theories. In particular, it is believed to be an effective theory 
for the deconfinement transition of SU(2) 30 ' 58 . The phenomenological viability of 
this approach will be fully addressed in another paper 52 . In any case, it is also 
commonly accepted that the computation of the effective potential for the Polyakov- 
loop requires the study of periodic scalar field theories, see e. g. 65 . Let us finally 
mention that direct studies of Yang-Mills theories using functional methods made 
recently major breakthrough, e. g. 69 > 34 , 57, 25, 27, 26,42,17^ w jvj snow below that our 
findings comfort these approaches. 

Let us now briefly introduce the challenges encountered when studying the bro- 
ken phase of the sine-Gordon model. While the UV of the broken phase is well under 
control being perturbatively renormalizable, the negativity of the /3-function for the 
amplitude of the model tells us that the IR will be strongly coupled, making the 
use of non-perturbative techniques such as functional renormalization mandatory 
(for reviews, see 59 > 10 > n ,56,28,22,60^ r^ ms f ea t ure has been studied using functional 
renormalization group in d = 2 51 , 49, 47, 44, 45,46,50,31 wnere a possible singularity re- 
lated to the convexity limit of the functional flows tends to appear in the IR. This 
same feature was also conjectured in d > 3 48 . The methodology of these studies 
consists in following the flow for the Fourier coefficients of the effective action, thus 
implicitly assuming that the relevance classification in the UV with respect to the 
Coleman fixed point still holds at least partially in the IR. We will detail why this 
assumption is incorrect in d = 2 due to the appareance of novel fixed points that 
can not be found using Fourier series, in agreement with previous findings 55 . The 
possible instability of the flow in d > 2 will also be shown to be a genuine one and 
thus of a different nature. 

While we believe that the non-pertubative behavior of the sine-Gordon model 
at LPA is now fully understood with the new results provided here, the inclusion 
of the wave function renormalization beyond the inclusion of a field-independent 
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anomalous dimension 46 at the non-perturbative level is still missing. This point is 
obviously of crucial importance in the case the sine-Gordon model as one expects a 
non-trivial anomalous dimension to yield the Kosterlitz-Thouless cross-over already 
at the perturbative level. The difficulty lies in the fact that the Effective Average 
Action (EEA) approach 66,67 which allows to define properly the next to leading 
order in the gradient expansion is based on the addition of a regulator that breaks 
explicitly the periodicity. As a result, the use of this RG scheme is non-trivial and 
some justifications have to be worked out in the case of the sine-Gordon model. 
On the other hand, the RG scheme of Wegner-Houghton (WH) 64 does not suffer 
from this plague as the fluctuations are computed explicitly with a loop expansion 
for some small parameter (see below). Its sharp cut-off nature does not allow the 
consistent definition of the wave-function renormalization but can provide at LPA 
reliable results that could be used as a reference for the EEA computations. 

1.1. Outline 

As a result from the previous considerations, we will adopt the following strategy. 

In section II, we introduce WH-RG and briefly discuss the main properties of its 
flow equation. We focus on the case d = 2 and point out the relevant quantities that 
characterize the flow as one expects the effective potential to be both convex and 
periodic at k=0 i. e. constant in both phases. We study perturbatively the location 
of the phase boundary and show the presence of a continuum of non-perturbative 
fixed points in the IR of the broken phase previously found in 55 . We then discuss the 
convergence property of the Fourier representation of the potential and show that 
only a global-non-Fourier-expanded numerical treatment can describe accurately 
the broken phase and prevent the appearance of singularities in the IR flow. The 
set of fixed points, which stabilizes the loop expansion and ensures a safe convex 
limit for k — > is then shown to be reached independently of the bare coupling u. 
We also extend the validity of the phase boundary to strongly coupled bare theories 
and demonstrate that the renormalized potential depends only on /3. 

In the section III, we derive the average action flow equation and project it 
on the LPA in d — 2. The evolution equation for the potential is then discussed 
and the impact of the regulator on the location of the perturbatively-establishcd 
phase boundary is discussed. In particular, we demonstrate that all the regulators 
provide the same phase boundary. The Fourier series convergence is also discussed 
and is shown to be even worse than in WH's case. It is then shown that for various 
regulators, all the features found in the previous section still hold. 

In section IV, we introduce the main features of the sine-Gordon model in d ^ 2. 
The UV perturbative scaling regime is exhibited and the phase diagram is discussed. 
The IR fixed points found in d = 2 are shown to diseapear and a fundamental 
instability of the flow in the IR for d > 2 is emphasized to be a genuine property. 
We then discuss the use of the average action and provide the analytical expressions 
that characterize the flow for various regulators. 
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Finally, we draw our conclusions in section V. 

To summarise, the present paper is adressing the following issues : 

• How the flow behaves in the deep IR of the broken phase and in partic- 
ular how the minimum selection is dynamically performed by the coarse 
graining. 

• How the perturbatively well-established phase boundary depends on the 
bare coupling u. 

• How the convexity of the effective potential is realized by the quantum 
fluctuations and the stability of the corresponding trajectories. 

• How universal are the trajectories in the deep IR. 

• How the addition of the regulator influences the previous results at LPA, as 
a preliminary step before the study of the wave function renormalization. 

• How the situation is changed in dimensions different than 2. 

2. Wegner-Houghton approach 

Let us start with a bare action defined at the UV cut-off k : 



where the field configurations 3>(x) have non-zero components only in the 
momentum-sphere of radius k. We are interested in the effective theory for the 
low-frequency modes <fi(x) which have momentum components only in the sphere of 
radius k — 5k once the fast modes (p(x) carrying modes in the shell [k — 8k, k] have 
been integrated out. Splitting the original field into its low and fast components 
= (f)(x) + ip(x) and performing a one-loop expansion for the fast mode ip yields 



where the Tr' holds for a trace in the momentum space on a shell of [k — 5k, k] and 
is thus proportional to 5k and ip is the saddle-point. The blocked action described 
by the previous flow equation can easily be proven to be the effective action, sec 
e.g. 5 . 

The presence of the saddle point (po m the flow equation opens the possibility 
to have tree-level renormalization 3,6 . This happens in particular in the IR of the 
broken phase for a </> 4 scalar field theory in the internal region between the two 
minima where tree-level effects provide the convexity of the effective potential taking 
into account inhomogeneous field configurations. Let us mention that even if this 
effect helps us to obtain a convex potential for k — > 0, this solution corresponds 
to a singularity of the flow equation with loop-evolution when the saddle-point is 
assumed to be trivial 2 ' 54 . Such a tree- level has been expected to occur in a similar 
fashion in the broken phase of sine-Gordon model in d=2 44 and in d > 3 48 in 
the concave regions using Fourier decomposition of the effective potential. We will 





(3) 
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show in the following that it is not the case for at least a part of the broken phase 
in d = 2, the whole IR trajectories being beyond our numerical capacities. 

Finally, let us mention that the flow equation (3) for the action is so far a one- 
loop approximation. At this point, one can notice that the higher loops terms that 
have been neglected carry at least two momentum integrals, each of them being 
proportional to 5k so that in the exact differential limit 5k — > 0, the flow equation 
(3) turns out to be exact. As the functional determinant is of little use for practical 
purpose and is in general unknown for an arbitrary background field cj>, one defines 
the Wegner-Houghton equation to be the projection of the evolution equation for the 
action on a constant background. As the blocking on sharp regions of the momentum 
space prevents us to define the next to leading order in the gradient expansion, see 
e. g. 15 - 59 7 the LPA is the only order of the gradient expansion that makes sense. 
The flow equation for the potential reads : 

where Kd is a geometrical factor given by : 

The equation (4) holds as long as the curvature of the action at LPA is strictly 
positive in order to maintain the validity of the loop-expansion. The limit of appli- 
cability/reliability of this equation thus reads : 

1 + V"{4>) - (6) 

where V" — V" /k 2 and will be referred to as the convexity limit of the blocked 
action in the following. 



2.1. UV perturbative regime 
2.1.1. Coleman Fixed point 

As we are looking for fixed points in the UV, let us define the flow equation for the 
dimensionless potential V in d — 2: 

k^ + 2V = -±-Lo 9 (l + V») (7) 

where we can make the confusion between the dimensionless and dimensionfull 
quantities <j> = <j> and /3 = ft due to a trivial scaling dimension. 

The Coleman fixed point is defined to be a perturbative UV fixed point for a 
periodic potential when no wave-function renormalization is taken into account. 
Taking the bare potential : 



V A = u A cos(f3(f>), m a /3 2 << 1, 



(8) 
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and linearizing the flow equation (7), one gets : 

dV 

47rfc— + V" + 8ttV = (9) 
ok 

where we neglected higher powers of (V") n = 0(u n ). This shows that the Coleman 
fixed point is characterized by the critical frequency : 



& = VStt (10) 

which reproduces the well-known result 20 . 

Defining f} r = j^-, the UV behavior in the weak coupling regime is easily estab- 



lished 



2 



I" 1 



2u fc (/3 r 2 - 1) (11) 



As a result, one expects that for /3 r > 1 the periodicity is unbroken due to vanishing 
of the dimcnsionless coupling with the flow while in the broken phase j3 r < 1 
the dimensionless coupling increases with the flow i. e. it exhibits UV asymptotic 
freedom. Taking into account that the flow generates during the blocking higher 
frequency modes, it is straightforward to check that the critical frequency for the 
higher Fourier modes n(3 is given by j3 r = 1/n within this linearization 44 so that 
when P r — 0.5, the Fourier mode of frequency 2/3 is marginal. 

On the contrary, in both phases, the dimensionfull coupling vanishes : 

= 2u k (3l > (12) 

This is the signal that one expects a dimensionfull effective potential to be a constant 
in both phases being the only function periodic and convex. We conclude that the 
dimensionfull quantities are not good quantities to characterize the theory. 



2.1.2. Order parameter 

As the scale is decreased in the broken phase, one expects at some point that 
our approximation of small coupling regime breaks down. The upper bound of the 
dimensionless coupling is given by u — » which corresponds to the appearance of 
soft modes with a single harmonic ansatz. Indeed, in such a case, the flow equation 
(7) is given by : 

k % + W = -^Log{l-p 2 ucos{P4>)) (13) 

where soft modes have appeared in the logarithm every ^p. In such a case, the 
linearization of the flow equation (9) is not valid anymore and one has to solve the 
full equation (7). 

The resolution for the Fourier expansion of the potential for the first 10 har- 
monics 46 and references therein showed that the flow in the broken phase seems to 
drive the potential to be Vk = —^<j> 2 repeated periodically in the concave regions. 
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This potential corresponds to the limit of applicability of the loop-expansion used 
to derive the Wegner-Houghton equation (4), making this scheme of renormaliza- 
tion questionable in such situation. Inspired from our previous considerations, one 
concludes that one should monitor the dimensionless curvature of the action at the 
middle of the concave regions where the convexity limit (6) has the best odds to be 
violated : 

S'^ n ) = l + V"^ n ), i = (14) 

At the perturbative level, this quantity has a /3-function of opposite sign as the 
one of the sine-Gordon amplitude each side of the critical surface i. e. it is UV 
asymptotically free in the symmetric phase and positive in the broken one : 

fc«(0„) = 2/3 2 £ fe (l-/3 2 ) (15) 

In addition, it provides a kind of effective order parameter as : 

• in the symmetric phase, this goes to 1 as the flow suppresses the coupling. 

• in the broken phase, this a priori goes to if the system goes unstable in 
the IR. 

This phenomenological order parameter has also the huge advantage to be easily 
computed in fRG, on the contrary to the topological susceptibility that has been 
demonstrated to be a natural order parameter for periodic theories 51 . If we are to 
make the connection with the SU(2) deconfinement transition, the curvature of the 
action has the opposite behavior than the one of the trace of the Polyakov loop : in 
the deconfined phase where the center symmetry is spontaneously broken, the trace 
of the Polyakov-loop is non-zero while it vanishes in the confined phase. 

2.2. IR of the broken phase 

2.2.1. Weakly coupled bare theory 

This study has been partially performed already in 55 but we reproduce here refined 
results for completeness. Let us first focus on the case where one starts with a 
perturbative initial condition /3 2 wa << 1, i- e. the UV scaling is given by equation 
(11) holds. Starting with a "weak" bare amplitude (3 2 ua = 0.01, the IR of the whole 
broken phase exhibits a continuum of fixed points for 0.90 > (i r > 0.55, see fig. 1. 
This line of fixed-points ensures the stabilization of the loop expansion and thus the 
safety of the use of the Wegner-Houghton RG. They correspond to an highly non- 
trivial recombination of all the Fourier harmonics and thus can not be seen using a 
Fourier truncation of the effective potential. This point will be further illustrated in 
section 2.2.2. A definite answer for lower frequencies than (3 r — 0.55 or larger than 
(3 r = 0.90 is unclear as one faces in the first case a very strong attraction of the 
convexity limit and a very low speed of the flow in the latter case. A study of the /3- 
function for the curvature illustrates well these difficulties, see fig. 2. We identify two 
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Fig. 1. The curvature at the middle of the concave regions (15) as a function of the scale k for 
/3 2 ua = 0.01 for increasing values of f3 r from 0.55 (bottom) to 0.90 (top) by steps of A/3 r = 0.01. 
The UV scaling is shortened for decreasing j3 r , in relation with cq. (11). 



different scaling regimes separated by a maximal speed of the flow, corresponding 
to the competition between the UV Coleman scaling and the IR fixed points. The 
perturbative UV scaling is shortened as /? is decreased, due to a larger speed of the 
flow of u, see (11) and (15). 

For f) r above 0.90, the transition between the two regimes is shallow and the 
/3-function is not enough suppressed relatively to the UV scaling to draw a clear 
picture. For f3 r below 0.55, the large bump in the /3-function signals the strong 
attraction of the singularity of (6) on the flow but the numerical errors prevent us to 
stabilize the IR scaling to establish undoubtedly -even if a strong suspicion is raised 
by the figure 1- the existence of an IR fixed point. An important aspect here is that 
the fixed points exhibited here are not only fixed points of the curvature S'^{4>) in 
<f> = <f> n , but also fixed points for the whole range of <f> € K. Let us illustrate it in the 
special case j3 r = 0.70, see figure 3. The change of sign of the /3-function around j| is 
easily understood when monitoring directly the curvature, as shown on the figure 4. 
It displays the curvature as a function of <f> for different scales. As it is well known, 
the spectrum of the broken phase of sine-Gordon is made of massive excitations, 
so that the dimensionless curvature in the convex regions is strongly enhanced as 
the scale is decreased below their typical mass scale. This point will be further 
discussed in section 5. In the same time, the curvature vanishes in the concave 
regions reflecting the progressive decoupling between the different homotopy classes 
of the dual XY model. These two effects in convex and concave regions are stronger 
as r is decreased. Finally, the IR dimensionless effective potential is plotted on 
fig. 5, where we see that the maxima of the curvature of fig. 4 corresponds to deep 
minima of the potential, providing an effective mass term, while the minima of the 
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Fig. 2. The /3-function of S"(<j> n ) in units of Afc as a function of the scale k for /3 2 Sa = 0.01 for 
increasing values of r from 0.55 (top in the UV, bottom in the IR) to 0.90 (bottom in the UV, 
top in the IR) by steps of A/3 r = 0.01. One can identify a power-law decay of the /3-function when 
reaching the neighborhood of the IR fixed points. 
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Fig. 3. The /3-function of the absolute value of the curvature \kdhS" (<j>)\ as a function of (j> for 
/3 2 ua = 0.01 for 0.7 . The vertical lines correspond to a change of sign of the /3-function : it is 
negative in the neiborhood of j|, positive elsewhere. For decreasing scales, the absolute value of 
the /3-function decreases uniformely for all <j>. 



curvature give us the upside-down parabola in the concave regions of the potential. 
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Fig. 4. The curvature of the action S" ((j>) as a function of <j> for four different scales for /3 r = 0.7 
in d=2. The peak (resp. flat) regions correspond to the convex (resp. concave) regions for the 
potential, see fig. 5. This shows that the peaks in the second derivative correspond to deep wells in 
the potential reflecting the progressive selection of a single vacuum, fluctuations around the latter 
being very massive. 
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-40 1 1 1 1 1 

Jt/p 2jt/p 3it/p 4ji/p 

Fig. 5. The rescaled fixed point potential /3^V with respect to <j> for for /3 r = 0.7 in d=2. The 
concave regions where the curvature is flat are given by the upside-down parabola — 10.998/3^ 2 , 
where 1 — 0.998 = 0.002 is the minima of the curvature at (f> = <j> n of fig. 4. 



2.2.2. Fourier expansion convergence 

Since the potential is periodic in the UV and since the equation (7) preserves pe- 
riodicity, one could have expected a Fourier series to provide a sensible expansion 
for the potential and try to limit ourself in following the scale evolution of the cou- 
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plings of the first Fourier components. Then, to ensure that we describe effectively 
the genuine solution, one could have checked the stability of such truncation when 
taking into account more harmonics. This strategy, followed so far in all the studies 
of sine-Gordon model is shown here to be misleading in the IR of the broken phase, 
in the same fashion as the polynomial expansion is in the broken phase of <f> -model 

2 

For the reasons mentioned previously, one solves the flow equation for the sec- 
ond derivative of the potential. In principle, as it is expected to be periodic and 
continuous, Dirichlet theorem should ensure us at least a simple convergence of the 
Fourier serie towards the genuine solution when increasing the number of modes 
taken into account. In the symmetric phase, this convergence is expected to be very 
fast as the lowest frequency is suppressed by the flow and thus has little possibility 
to generate higher frequency modes which are themselves even more irrelevant. This 
is indeed what is found in 44 and even a single harmonic ansatz for the potential 
could do the job. 

On the contrary, the broken phase is more involved. First, the original lowest 
frequency is relevant and is enhanced by the flow, generating higher frequencies 
which may be themselves enhanced by the flow (if f3 r < 0.5) but are in any case less 
irrelevant as in the symmetric case. In addition, the higher frequency modes have 
an enhanced weight in the second derivative of the potential, their amplitude being 
the one of the potential times p 2 for the p-th mode. 

But the largest source of error and the origin of the bad convergence of the Fourier 
expansion comes from a kind of Gibbs phenomenon in the concave regions of the 
potential. To illustrate this, let us define the Fourier expansion of our numerical 
solution : 

n 

V"{4>)=Y J C P cos{p^) (16) 
P =i 

where the c p 's can be reconstructed from the usual identity : 

c P (k) = f - f + ' 3 14» cos{pp<t>)d<t> (17) 

Let us illustrate for the special case r — 0.7, (3 2 ua = 0.01. This initial condition 
corresponds to a trajectory entering the broken phase. Starting from the UV, the 
curvature weakly decreases and is then strongly attracted by the convexity limit (6) 
in the IR. For the second derivative of the potential, it amounts in the appearance 
in the IR of flat regions of constant value ~ —0.998 where the potential is concave 
and peaks in the convex regions, see figures 4-5. Let us notice that this shape can be 
obtained only by the infinite resummation of all Fourier components as the higher 
order harmonics are IR enhanced. On figure 6, one sees that when the original 
coupling becomes large enough, it triggers the build up of the higher modes, one 
after another. Since the sign of their contribution depends on the parity of the 
modes, the fixed point solution is a finely-tuned competition of all the modes. This 
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Fig. 6. The absolute value of the reconstructed Fourier coefficients of the second derivative of 
the potential defined by eq. (17) as a function of the scale k for /3 r = 0.7 and «a = 0.01. The even 
modes are negative and the odd ones positive. 

situation is comparable to the broken </> 4 theory where we know one needs the exact 
resolution for the effective potential without any polynomial ansatz to recover the 
Maxwell-construction when projecting on the polynomial couplings 54 . 

Even if the transition between the convex and concave regions has to be con- 
tinuous, it is in practice rather sharp due to the logarithmic enhancement, cf (7). 
As a result, the Fourier serie tries to mimic a change in the second derivative that 
resembles a step- function, see fig. 4 (note the log y-axis). The serie of the cosine- 
polynomial that are a set of continuous functions can not hope to converge uni- 
formly towards an almost discontinuous function near the frontiers between convex 
and concave domains of the potential. Thus the partial sums of eq. (16) overshoot 
the second derivative of the potential and this overshoot does not die out 70 . As a 
result, the truncation of the flow equation on the lowest frequency modes introduces 
spurious oscillations and is not able to stabilize before the convexity limit (6). This 
is shown on figure 7 where the Fourier oscillations are shown to be violating the 
degeneracy criterion (6) while the genuine solution does not. Let us mention that 
this overshoot is unseen using the Parseval identity : 

h = 7r l^"WI 2 # - E M fc )l 2 = lim 5 « ( 18 ) 

a P =i 

see figure 8. The discrepancy is shown in the latter figure to die out as the integral 
of the genuine solution is not related to its continuity properties. 

Due to this phenomenon, the Fourier expansion in the broken phase is of dan- 
gerous use, at least in the IR. From our point of view, this explains why the other 
studies 46 and references therein incorrectly predicted the appearance of the singu- 
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Fig. 7. The curvature of the action S" as a function of <j> for the scale for k = 0.79 ■ 10~ 4 at 
which the fixed point solution has been reached in d=2 close to the edge of the fiat region. The 
Fouricr-Dirichlet sum stands for the partial sum of the Fourier serie up to order n = 17, see eq. 
(16). The genuine solution is the solution coming from the resolution of the flow equation without 
Fourier expansion. 




Fig. 8. The relative error of the Parseval sums I Sn Ik Ik I as a function of the scale k, see eqn. (18). 
The first harmonic is a good description only in the UV perturbative regime, up to k ~ 0.1 and 
represents less than 40% of the genuine solution in the IR, when the fixed point is reached. The 
relative error decreases monotonously as n is increased. 

larity and failed to see the fixed points line as the overshoot drove the flow in the 
unstable region for the flow equation. The absence of the instability -the breakdown 
of loop-expansion- is of crucial importance in our opinion. If it was a genuine in- 
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stability, it would mean first that the RG scheme used to find it was inapplicable 
in this case. But this would also signal that the related path integral sees equal 
contributions of all field-configurations in the concave regions, a situation that can 
not be handled analytically nor numerically so far. Having clarified this point, we 
present in the following the solutions of the flow equation without using any Fourier 
ansatz, see appendix Appendix A. 

2.3. Universality in the broken phase 

We showed in the previous section that at least a part of the broken phase was having 
IR non-perturbative fixed points when the bare theory is chosen to be perturbative. 
A natural question is thus to determine in which extent the coupling ua has to 
be small and study the influence of its bare value on the existence of the IR fixed 
points and -if still exist- how the way they are reached changes. 

It is well known when the wave-function renormalization is taken into account 
perturbatively that for large enought couplings, the presence of a separatrix in the 
phase diagram places the system in a cross-over regime 33 and even simple ansatz 
in fRG reproduces it well, 46 . 

On the other hand, if the initial coupling is really large, the perturbative lineariza- 
tion of the flow (9) breaks down as non-perturbative effects arise at the cut-off scale 
and no simple analytical method can locate the phase boundary. In such situation, 
an LPA study might teach us something valuable as we have in any case no argu- 
ment to put forward to claim that we are in the cross-over regime that depends 
crucially on the wave-function renormalization. Starting with /3 2 ma < 1 i.e. close to 
the break down of the loop-expansion in the UV, one expects the strong enhance- 
ment in the UV of the higher frequency modes in the same way as in the IR of the 
broken phase. As these modes have a relative weight larger in the second derivative 
of the potential -the p-th mode is p 2 enhanced- they can compete in the /3-function 
of the curvature with the original mode. Since this happens already in the UV, one 
should expect non-trivial effects to appear and even possibly a change in the phase 
diagram. Let us show that it is not the case at LPA order. As usual, we illustrate 
the results with the case j3 r = 0.7. On figure 9, one sees that in the UV the flow of 
the curvature tends to restore the symmetry for large enough bare coupling, at odd 
of the perturbative relation (15), then the curvature saturates and finally reaches 
an IR fixed point. For smaller values of the bare coupling, the /3-function remains 
positive all along the trajectory in agreement with the perturbative expectations, 
eq. (15). The power-law scaling found in the IR is a universal quantity, characterized 
by an exponent independent of u\ (but dependent on (3 r ). 

More interestingly, the figure 10 shows that not only the approach to the fixed 
points is universal, but also the value of the curvature itself is independent on the 
bare coupling. Let us mention that this generalized universality spreads for the whole 
range of </>, as exhibited on figure 11, where the relative error on the IR curvature 
for different values of ua is shown to be very small. The specific shape of the plot is 
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Fig. 9. The /3-functions for the curvature S"(<p n ) in units of Ak for fi r = 0.7 for /3 2 «a £ [0.1 
0.99]. 



driven by the shape of the not-exactly zero /3-function of the curvature at finite k. A 




-itr 4 io" 3 icr 2 io" 1 0.9 0.92 0.94 0.96 0.98 
Fig. 10. The curvature S"{4>„) with respect to the scale k for /3 r = 0.7 for /3 2 u A G [0.1 : 0.99]. 



study of the Fourier couplings demonstrates immediately that the UV enhancement 
of the curvature for strong enough bare coupling is due to the dominance of the 



February 8, 2012 1:37 WSPC/INSTRUCTION FILE sinegordonl 



16 V. Pangon 



6.0 10" 5 




Fig. 11. The relative error on the deep IR curvature S" '(<f>) with respect to <j> for f} r = 0.7 
for /3 2 ua = 0.8 and 0.9. The reference is chosen to be the fixed point curvature computed for 
P 2 u A = 0.01 



higher frequency modes that are strongly generated in the UV, see figure 12. The 




10" 4 10" 3 10" 2 10~ 1 



Fig. 12. The Fourier coefficients c p of definition (17) for /3 r = 0.7 for /3 2 u A = 0.01 and 0.99. The 
fixed point structure is the same in both cases, but the flow history is different. In the strong bare 
coupling case, the higher frequency modes arc generated already in the UV. 

property of generalized universality -the presence of the same IR fixed points at 
strong and weak coupling- has been checked to occur for any j3. r > 0.55. For lower 
r , the problem of facing the too strong attraction of the instability (6) is obviously 
still present. Noticing that the final effect of choosing a strong bare coupling amounts 
only in reaching faster the same IR fixed point, one can reasonably conjecture that 
for 1 > (3 r > 0.90 the fixed point structure is preserved and allow us to find the 
fixed point value for the curvature for 0.985 > (3 r > 0.90 using a strongly coupled 
bare theory. As a result, one can summarize the set of fixed points of the broken 
phase in figure 13. A clear answer about the sectors (i r < 0.55 and (3 r > 0.985 is 
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beyond our numerical capacities. 




0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 



Fig. 13. The fixed point structure of the broken phase i. e. the IR value of the curvature as a 
function of r . The values for r > 0.90 come from the conjecture that strongly coupled bare 
theory also converge to the same IR fixed as the weakly coupled one. 



3. Average action 

As one knows that the effect of the renormalization of the temperature of the dual 
XY model is strong, one needs to capture the main features of the sine-Gordon 
model to take into account the wave-function renormalization that amounts here in 
the effective change with the scale of the /3-frequency. Unfortunately, the Wegner- 
Houghton RG which belongs to the class of RG with sharp regulator does not make 
sense at the next to LPA order. To overcome this difficulty, one needs another 
RG scheme with a smooth regularization. The average action method 67 , being an 
"exact" scheme i.e. not based on any kind of approximation is the best candidate. 
We study in this section the price to pay for the use of the regulator at LPA where 
we can compare it to Wegner-Houghton as a preliminary work before addressing 
wave-function renormalization 

Let us start with the bare action 5a [<f>]. The goal is to build a functional that 
interpolates between the bare action and the full regular effective action by a pro- 
gressive turning on of the fluctuations with the help of a regulator playing the role 
of a momentum dependent mass term Mfc[</>] = \ J (j)Rk4>- The partition function of 
this modified theory now reads : 

/S A + M k +J4> 
V4>e- *" (19) 

The generating functional of the connected n-points functions is thus : 

W k = -hLogZ k (20) 
Let us define the so-called average action to be : 

T fc = minj {W k - J(j> c i} - M k [<j> c i] (21) 
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This shows that T+M is a convex quantity, being the result of a Legendre transform. 
One can then derive the flow equation for the average action by formal functional 
manipulations 67 where only the regularity of the Legendre transform (21) is needed. 

(2) 

This assumption amounts in having T) ' + R k > for any scale k during the flow. 
The resulting flow equation is then given by : 

, R, \ 

(22) 

Let us now turn our attention to the sine-Gordon model particular case, where 
Sa[4>] is a periodic functional in the field <j). The n-points functions computed from 
Wk or Zk a priori correspond to a theory where periodicity is explicitly broken by the 
addition of M k . On the other hand, one can gain some insights on these expressions 
using a one-loop evaluation of the partition function Zk of eq. (19) which yields : 

T fe = SM + ^TrLogiS'M+R^ + oih) (23) 

where the regulator contribution at tree-level has cancelled out. This loop expansion 
is certainly valid in the UV as the regulator dominates and behaves as a large mass 
term 67 . This shows that the initial condition for (22) can be safely chosen to be 
periodic. Then the flow equation preserves the periodicity as the regulator is involved 
only in r^ 2 ) + Rk i. e. it adds an effective mass to the fluctuations. This situation, 
comparable to the massive sine-Gordon model 45 , drives us to conclude that no 
further non-periodicity than the one from the regulator (that will eventually vanish 
at k — > 0) will appear. In particular, this allows us to study the second derivative 
of the average action with respect to as a periodic function of <f>. 



3.1. UV perturbative regime 

The Local Potential Approximation of the equation (22) reads in d dimensions : 

ov k i r d d kR 

dk ~ 2 J dP R + pi + V" 
It is convenient to express the regulator as : 

R k =p*r(^), kd k R = -2p'(^) (25) 



where the function r is constrained by the regulator prescriptions 



67 



lim r(y) = oo, lim r(y) = (26) 

y— >0 y— ¥oo 



This allows to rewrite (24) in d = 2 using the dimensionless variable y = £ 

s y 2 r'{y) 

yr(y) +y + V" 
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Following the same strategy as in the Wegner-Houghton case, we can linearize the 
flow equation in the UV : 



4?r J q r(y) + 1 

+ V"±- [°°dy , t~~t— — ~~ttj (28) 
47ri y (r(y) + l) 2 V ; 

where the constant field independent contribution of the right hand side of the first 
line is a left-over of the regulator and can be discarded here as a shift of the vacuum 
energy. From the linearized flow in the UV, eq (28) one can check that the Coleman 
frequency is the usual one: 

fi= ( * r^w r'.g. (29) 



8^ J (r(y) + l) 2 

This result has been obtained independently of the choice of the regulator and its 
parameterization, only using the general properties eq. (26). The considerations 
from the WH case generalize immediately and the harmonic cos(/3<f>) is found to be 
irrelevant for /? > /3 C , and relevant for /3 < /3 C . 

Let us now introduce the regulators that will be used in following. The regulator 
(sometimes called Litim's 36 ) is a very popular choice as the loop-integral (27) can 
then be computed analytically : 

n(y) = Q-i)e(i-i/) 



kd k V + 2V=- I —^ 7ri (30) 



2 1 

If we are to prepare the inclusion of wave-function renormalization, we will also need 
a smooth regulator, the Litim regulator having trouble with gradient expansion 19 . 
We pick for simplicity the power-law defined by : 

r P (y) = a p y~ h ?, (a p , b p ) E M+* x [1 : oo] (31) 

where the loop-integal (27) can be analytically performed, see 50 and appendix Ap- 
pendix C. This regulator for the particular choice b p = 1 i. e. Rk — k 2 is actually 
a functional Callan-Symanzik RG scheme 38 ' 7 ' 4 and the corresponding flow equa- 
tion can be recasted into the Wegner-Houghton equation (7), see 50 and appendix 
Appendix B. This alternative derivation for Wegner-Houghton from a scheme using 
a regularization from massive fluctuations provides us even more confidence that 
average action is able of dealing with periodic theories. 

To summarize, we showed that even if the average action scheme introduces a 
mass term to regulate the fluctuations, this insertion is actually soft as this does not 
break explicitly the periodicity of the initial condition and the massive fluctuations 
are periodic. The phase boundary established in the UV was shown to be preserved 
by the addition of any regulator and we can hope that the IR is mostly untouched 
due to the vanishing of the regulator. 



February 8, 2012 1:37 WSPC/INSTRUCTION FILE sine'gordonT 



20 V. Pangon 

3.2. IR of the broken phase 

3.2.1. Weakly coupled bare theory 

Let us redo the d = 2 study performed in section 2.2.1 for Wegner-Houghton RG 
with Litim's regulator for simplicity. We will see after that the qualitative and most 
of the quantitative results will spread to other regulators. Starting with a weakly 
coupled bare theory (3 2 ua = 0.01, we scan again the phase diagram for various f3 r . 
The equivalent quantity of the dimensionless curvature for Litim regulator is : 

S"(<f> n ) = 1 + V"(<j> n ) (32) 

which tests this time the possible singularity of the Legendre transform (21). This 
will be further discussed in the next section. The figure 14 exhibits the general ten- 

1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


10 8 10 7 10 6 10~ 5 10~ 4 10~ 3 10 2 10 1 10° 

Fig. 14. The curvature S"{<f> n ) with respect to the scale k for /3 2 «a = 0.01 for increasing values 
of p r from 0.40 (bottom) to 0.90 (top) by steps of A/3 r = 0.01. 

dencies of the flow for various j3 r . The first obvious remark is that the curvature still 
exhibits a continuum of fixed points. This feature is well preserved by the addition 
of the regulator as it could have been expected, being a phenomenon occurring in 
the IR where the regulator vanishes. The second remark is that the fixed point value 
of the curvature is larger as in Wegner-Houghton RG (note the linear scale for the 
y-axis). The reason is that the regulator which behaves as a mass term suppresses 
a part of the fluctuations that were shown to decrease the curvature in Wegner- 
Houghton picture. Finally, one can notice that we are able to stabilize a fixed point 
for lower f3 r values, namely around (3 r — 0.4 (one can stabilize few lower trajectories 
not displayed for clarity). Realizing that the curvature is still "large", we discover 
here that the convexity limit 1 + V" = is not the limiting phenomenon as it 
was the case in Wegner-Houghton picture. This will be shown to be related to the 
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non-differentiability limit of the running potential in the convex regions when using 
Litim regulator. These considerations are further supported by the /3-function of the 




Fig. 15. The /3-function of the curvature S"(<j> n ) in units of Ak with respect to the scale k for 
/3 2 ua = 0.01 for increasing values of f3 r from 0.40 (top in the UV, bottom in the IR) to 0.90 (top 
in the IR, bottom in the UV) by steps of A/3 r = 0.01. 

curvature, see figure 15. The general behavior of the /3-function is the same as the 
one obtained in Wegncr-Houghton namely the existence of a maximal speed deep- 
enough in the broken phase corresponding to the cross-over between Coleman and 
IR fixed points. The IR scaling regime when entering the fixed point neighborhood 
is also given by a power-law decay. Interestingly, the power is different to the one 
obtained in Wegner-Houghton case. The figure 16 shows that the UV scaling regime 
is the same, being given by perturbation theory in both cases. Around k = 10~ 2 the 
flow is entering the cross-over regime and Litim and Wegner-Houghton flows start 
to diverge. The latter is strongly attracted by the convexity limit while the former 
still feels the influence of the regulator suppressing the fluctuations. As a result, the 
trajectories end at different fixed point values and the fact that the power of the 
decay towards the fixed point are slightly different shows that the regulator impact 
at higher scales is not totally washed out. 

3.2.2. Fourier expansion convergence and convexity issue 

We saw in the Wegner-Houghton case that the Fourier serie had a bad convergence 
in the IR of the broken phase due to the appearance of large flat regions of the 
second derivative of the potential corresponding to a dynamical Maxwell-cut in 
the concave regions. The loop-expansion has been shown to be stabilized before 
the occurrence of the instability for a part of the phase diagram /3 r > 0.5. Using 
average action, the IR of a scalar broken phase in the concave region might also face 
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Fig. 16. The curvature and its /3-function in units of 10 4 Ak with respect to the scale k for 
/3 2 u A = 0.01 and /3 r = 0.70 

an instability. This instability is not related to the degeneracy of the average action 
as the derivation of the flow equation does not rely on a loop-expansion but comes 
from the Legendre transformation (21) which needs T y k '+Rk to be invertible i.e. 
r^. 2) + Rk > 0. At LPA, the convexity limit is thus given by 62 : 

fc 2 + V" = (33) 

where k — min p {p 2 + Rk(p)}- This limit has been found to be strongly IR at- 
tractive 62 in the internal region of a broken </> 4 model in d = 4 giving rise to 
inhomogeneous tree-level contributions below the scale of appearance of this would- 
be instability. In such a case, the IR potential in the concave regions is given by the 
upside down parabola V k = — \§ 2 and thus depends explicitly on the regulator. 
This raises the question of the possibility of removing the regulator for k — > in the 
internal region. Intuitively, one realizes this problem noticing that the regulator is 
a dimensionfull quantity so that it needs the appearance of another mass scale to 
decouple and then vanish. In the internal region, there is no other physical scale so 
that this decoupling can not occur at any finite scale k. The regulator being the only 
available scale in the non-convex region drives the flow and is itself considered as the 
typical momentum scale of the inhomogeneous configurations dominating the IR of 
the broken phase. Another way to say this is to notice that while the dimensionfull 
potential goes to in the internal region -reproducing Maxwell construction-, the 
dimensionless one reaches a fixed point V = — ^-<j) 2 corresponding at the convexity 
limit where k = k 2 k 2 . From these considerations, it is now obvious why the average 
action is able to provide the jump in the susceptibility at the border of the concave 
region : in the internal region the flow continues to stick to the regulator while on 
the contrary, the regulator has decoupled due to the presence of a mass term near 
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the minima (the Higgs mass) . The picture is even stronger in dimensionless quanti- 
ties, when the second derivative reaches the fixed point — k 2 in the internal region 

2 

and is quadratically enhanced in the external one (~ -£#■)■ 

Let us now turn our attention to the sine-Gordon model. From similar consid- 
erations than the ones used in WH case, the well suited order parameter is : 

S"(<f> n ) = ~k 2 + V"(^ n ) (34) 

where <ft n in defined by eq. (14). Due to this jump in the susceptibility, the Gibbs 
phenomenon invoked in WH case is even more predominant being typically e.g. a 
power- law enhancement using Litim regulator (cf. eq. (30)) to be compared the the 
logarithmic one in WH (cf. eq. (7)). So that we expect the Fourier serie convergence 
to be even worse in the average action case. Let us illustrate this for j3 r = 0.7 in 
d = 2 using Litim regulator for which k = 1. In Figure 17, one sees that in the IR, 

10 3 
10 2 
10 1 
10° 

10- 1 
io- 2 
io" 3 

Jt/p 27t/p 3n/p 4n/p 

Fig. 17. The curvature S"(<f>) with respect to <j> for f3 r = 0.7 in d=2 (continuous line) for k = 10~ 4 . 
The peak (resp. fiat) regions correspond to the convex (resp. concave) regions for the potential. 
The dashed line corresponds to the rcscalcd potential 0^V at the same scale and shows that the 
peaks in the second derivative correspond to deep wells in the potential reflecting the huge selection 
of a single vacuum. The concave regions where the curvature is flat are given by the upside-down 
parabola — |O.88/3^0 2 , in relation with the fixed point value of fig. 16. 

the wells in the potential selecting one minimum among the others is much stronger 
than in WH (note the logarithmic scale). As discussed previously, this is given by 
the jump of the susceptibility provided by the regulator. As a result, the comparison 
between genuine result and the Dirichlet partial sum provided by figure 18 shows 
a much larger discrepancy compared to WH case. Obviously, the deeper we are in 
the broken phase, the stronger is the selection of the minima. This consideration 
shows that the limiting (3 r where we are able to find a fixed point is not given by the 
convexity limit but rather by the limit where in the IR the selection is so strong that 
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Fig. 18. The curvature of the action k 2 + V" with respect to <j> for the scale for k = 0.25 ■ 10~ 3 at 
which the fixed point solution has been reached. The Fourier-Dirichlet sum stands for the partial 
sum of the Fourier serie up to order n = 15, cf. eq. (16). The genuine solution is the solution 
coming from the resolution of the flow equation without Fourier expansion. 



it amounts in a numerically non-differentiable interaction. To illustrate this point, 
let us study the /3-function for the whole curvature of figure 19. While the flow is 
clearly suppressed in the whole concave regions as the figure 15 was suggesting it, 
the /3-function is still strong in the convex regions when the flow is stopped by the 
non-differentiability of the interaction. On the other hand, we can notice that the 
/3-function is decreasing fastly in the IR so that there is good hope that it will also 
reach a maximal value. We will see in the next section that the fact that this small 
scale dependence at the end of our flow has no impact on our conclusions. 

3.3. Universality in the broken phase 

Having found that the average action allows us to find the same set of fixed points 
as in the Wegner-Houghton case, it is tempting to test if the universality found 
previously still holds. This is actually far from obvious as picking a strongly cou- 
pled bare theory changes its relative weight to the regulator in the UV. We will 
demonstrate here that this difficulty does not occur. 

The figure 20 shows the /3-function of the curvature of the theories for f3 r = 0.7 
and various ma- For couplings large enough, the /3-function of the curvature also 
changes sign along the trajectory and exhibits a common power-law decay towards 
the fixed point independently of u\. The figure 21 shows that the fixed points values 
are also independent of u\, in the same fashion as in WH case. The figure 22 allows 
us to check that the whole curvature, not only at the middle of the concave regions, 
is the same for different large bare couplings. 

Obviously, all these results can be generalized to other f3 r values and as a strong 
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Fig. 19. The /3- function of the absolute value of the curvature | kd^ S" (</>) | with respect to <j> at 
the end of the flow for /3 2 ua = 0.01 for r = 0.7 for two different scales. A log-scale is used for 
the x-axis and the value <j> = |y is meant to be reached only asymptotically. The change of sign of 

the /3-function in ^ and is hidden here for clarity as it is reduced to a vertical line. 




bare coupling amounts only in a faster convergence to the fixed point, one can 
establish its presence for 0.90 < r < 0.98. 

Finally, one can summarize the fixed point structure of the broken phase in 
figure 23 where the difference with Wegner-Houghton is clear. Let us notice that 



February 8, 2012 1:37 WSPC/INSTRUCTION FILE sinegordonl 

26 V. Pangon 



u A =aio "- — h — — ' ""' 

u.=o.2o —x— y 

u A =0.30 x /> 
U A =0.40 b // 
U A =0.50 -— - / // 

u A =0.70 - - / ' / b 

u A =0.90 ----- / x*ifc 
u A =0.99 / ; /Ms, 


i 

«xxx* x-x x-x-x-x-x x-xx x-x-x x-x-x-x-x-x i 
Kxxx-x^xxxxxxxx-x-x-x-x-x-x-x-x-x-xxi 


3 BOB □□BQBOaaBBBDBBOaD B-B B O-E 
■■■HHI 1 1 1 IHl | ■ ■■ M Hfl , { 

" 


1 









0.9 0.92 0.94 0.96 0.98 



Fig. 21. The curvature S" (</>„) as a function of the scale fc for r = 0.7 for /3 2 u A e [0.1 : 0.99]. 
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Fig. 22. The error on the fixed point curvature for f} r = 0.7 for ua = 0.80 and 0.90 in d = 2 
relatively to the curvature computed for /3 2 ua = 0.01 at the end of the flow. A log-scale is used 
for the x-axis and the value (j> = ^ is meant to be reached only asymptotically. 



the fixed point curvature as a function of f3 r has a clearly different behavior in the 
neighborhood of r ~ 1. 
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Fig. 23. The fixed point curvature using average action with Litim regulator in function of f) r . 
The Wegner-Houghton results are displayed for comparison. 

3.4. Influence of the regulator 

One of the final goal of this study being to prepare the situation for the inclusion of 
the wave-function renormalization, we have to be sure that the results established 
previously still hold for a smooth regulator and are not an artifact of the regulator 
choice. In this sense, this section can be seen as the first steps of an optimization 
study. 

As expected the flow equation to solve is much more time consuming being an 
intcgro-partial differential equation as the momentum integral (24,27) can not be 
performed analytically in general. Let us present first the results for a power-law 
regulator given by (31) for a p = 1 and b p — 2, see 50 and appendix Appendix C. In 
such a case, one has to use k = 2 in the definition of S"(<frn), see (34). 

The universality found using Litim regulator is still present, namely the inde- 
pendence on the bare sine-Gordon amplitude of the IR /3-functions, see fig. 24. The 
independence of the bare sine-Gordon amplitude is easily established as the curva- 
ture trajectory ends at the same value for the both weak and strong bare coupling, 
see figure 25. 

As we saw previously, the whole curvature sees a set of fixed point that are the 
product of a finely tuned balance between the convex regions where the curvature 
is strongly enhanced and the concave ones where the curvature being attracted by 
the convexity limit tends to vanish, see figure 26. At this stage, one can notice that 
the curvature is indeed the good order parameter as in this case k 2 = 2 but the 
curvature has still the same order of magnitude of values as in Litim's case. Obvi- 
ously, the precise value depends on how the regulator suppresses the fluctuations. 
An obvious preoccupation is to study the dependence of the fixed-point solutions 
when changing the regulator and in which cases we can get rid of this spurious 
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Fig. 24. The /3-function of the curvature at the edge of the concave regions with respect to the 
scale k near the fixed point for /3 r = 0.7 and /3 r = 0.8 for /3 2 Sa = 0.01 and /3 2 «a = 0.80. This was 
computed with b p = 2. 




Fig. 25. The curvature S"(4> n ) with respect to the scale k near the fixed point for /3 r = 0.7 and 
I3 r = 0.8 for /3 2 w A = 0.01 and /3 2 u A = 0.80. This was computed with b p = 2. 



effect. Using parameterization (31), one finds : 



~k 2 = (flp^p-l)) 1 



1 + 



b P -I 



This converges to 1 in the large b p limit slowly: 

£2^1 , l °9(b P ) 



(35) 



(36) 
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Fig. 26. From top to bottom : the fixed point curvature S"(<f>) for fi r = 0.7 for /3 2 ua = 0.01 and 
the corresponding relative error times 10 3 on the curvature computed for /3 2 «a = 0.80. The last 
plot is the /3-function for the whole curvature \kdkS" (rf>)\ times 10 s for /3 2 ujy = 0.01. This was 
computed with b p = 2. 

and diverges in the large a p limit which amounts in rejecting the convexity singu- 
larity at k = 0. We are obviously not interested in the latter limit as the IR fixed 
points arc found at finite distance of this singularity. The figure 27 shows that the 
discrepancy between the fixed-point solutions for different power b p grows as we get 
closer to the Coleman fixed point. For the related /3-functions, the behavior is the 
oppposite, as expected, see figure 28. 
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Finally, let us compare the shape of the fixed point curvature for the different 
RG schemes and regulators presented here, see figure 29. In the convex regions, the 
curvature is strongly enhanced in all cases. While in Litim's case this amounts in an 
interaction that resembles a non-differentiable one, one can smoothly match the cur- 
vature given by average action with the one given by Wegner- Houghton when using 
the power-law regulator for increasing powers. This amounts in concluding that the 
would-be non-differentiability is not a genuine one, but rather an artifact of the non- 
analyticity of Litim regulator. In the concave regions, the difficulty of handling a 
system close to the limit of applicability of the loop-expansion for Wegner-Houghton 
is soften by the use of average action as its distance to the would-be singularity of 
the Legendre transform lies much further. 

4. Other dimensions 

Let us now focus on the case where the dimension is not 2. As we saw already, 
the dimcnsionfull flow equation preserves the periodicity in any dimension for all 
RG schemes. The dimcnsionfull effective potential can thus be Fourier-expanded at 
least in the UV perturbative regime : 

oo 

V k (cl>) = J2un(k)cos(nf3ct>) (37) 

n=l 

where the dimcnsionfull frequency (3 is a constant of the flow. On the other hand, 
due to the finite scaling dimension of the field, the dimensionless flow equation 
breaks periodicity. This amounts in Wegner-Houghton RG to : 

k%+dV- d -^W' = -^Log(l + V") (38) 
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Fig. 29. The fixed-point curvature at the end of the flow with respect to <fi for 
/3 2 «a = 0.01 in d = 2 for different regulators and Wegncr-Houghton for comparison. 



0.7 for 



where the dimensionless effective potential can then be written in terms of dimen- 
sionless variables: 



y(0) = ^u„(fc) COS (n/3(fc)0) 



(39) 



n=l 



The dimensionless frequency is now scale dependent 48 : 

/3 2 (fc) = p 2 k d - 2 (40) 

The first obvious conclusion is that the Coleman frequency does not exist in d ^ 2 
as the dimensionless potential is not periodic : even if the amplitudes u n were seeing 
a fixed-point, the frequency would still run. As a result, any perturbative periodic 
bare theory is driven by the the gaussian fixed point. Let us linearize the flow 
equation in the UV weak coupling limit : 



(41) 



Thus the UV perturbative flow is driven by : 

kd k Un = u n d (n 2 p 2 r k d - 2 - 1) 

u n (k) = u n (A) (-) e «" 2 UaJ - 1 ) 

where we defined in a similar fashion as in d — 2 the reduced frequency (3 2 — ^jfl 2 . 
By abuse of notation we will note : 



(42) 




February 8, 2012 1:37 WSPC/INSTRUCTION FILE sine'gordonT 



32 V. Pangon 

as it reproduces the Coleman frequency in d — 2 and separates two different UV 
behavior of the /3-functions in d ^ 2. In the following, we use the UV cut-off A = 1 
for simplicity as it provides (3(A) = (3. 



4.1. No symmetric phase in d > 2 

Let us first study the situation for d > 2. For (3 r < 1/n, the /3-function of the mode 
n is always negative and exhibits UV asymptotic freedom. This corresponds to the 
beginning of a broken phase scaling, where the coupling grows with the flow. At some 
point, the couplings are large enough to break the perturbative assumption needed 
to linearize the flow and the full flow equation has to be solved non-perturbatively 
to conclude about the continuation of the flow. 



10° 
10 6 
10 4 
10 2 
10° 
10' 2 
10' 4 

io- 6 

10" 8 



S" p.=0.90 
S"P r =1.10 
S" p r =2.00 
u, p.=0.90 
u, B,=1.10 
i r =2.00 




Fig. 30. The flow of the curvature S"(<j> n ) and the first Fourier harmonic «i in units of /3 2 with 
respect to the scale k for different /9 r with /3 2 Si(A) = 0.001 in d = 4. On the right hand side, the 
curvature value is given by f» 1 due to the very weak bare coupling for the three different values 
of p r . 



For j3 r < 1/n, the /3-function of the mode n starts with a positive value and then 

i 

changes sign at the scale k n — {^rj^j ■ In such a case, the trajectory drives the 
system closer and closer to the Gaussian fixed point as the frequency is increased 
to be finally reppeled to the broken phase in the IR. The UV flow of the curvature 
is then given only by the first harmonic as the higher ones, generated by the small 
non-linearites of the flow equation, are relatively strongly suppressed. These consid- 
erations are summarized in figure 30. As a result, one can conclude that at LPA, the 
bare single harmonic sine-Gordon model is always driven by the perturbative flow 
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to its broken phase. At the non-perturbative level, the IR behaviour of the figure 30 
tells us that the broken phase behavior is not stabilized anymore by non-vanishing 
curvature fixed points as in d — 2 and thus yields the appearance of the instability 
i. e. a degeneracy of the curvature signaling the breakdown of the loop-expansion. 
This instability, predicted to occur in 48 using Fourier series, is demonstrated here 
to be a genuine one, being established with a global algorithm. 

4.2. No broken phase in d < 2 

The previous discussion still holds in the same fashion as before except that having 
d < 2 reverses the behavior. More specifically, for f3 r < 1/n, the n-th sine-Gordon 
harmonic starts growing up to k = k n and then decreases. As a result, the per- 
turbative scaling regime is shortened as the higher modes are enhanced and pile 
up during the evolution. Provided that (3 r is not too small, one can nevertheless 
preserve a perturbative trajectory that will at the end of the flow suppress all the 
couplings. For f3 r < 1/n, the flow monotonously suppresses the couplings. These 
considerations are illustrated on figure 31 where we show that even when the flow 
reaches the non-perturbative neighborhood of the instability, it finally goes toward 
the gaussian fixed point. 




Fig. 31. The flow of the opposite of the potential curvature —V" ((fin) = 1 — S" (<fi n ) and the first 
Fourier harmonic ui in units of /3 2 with respect to the scale for different /3 r with /3 2 «i(A) = 0.001 
in d = 1. We see that starting with a small f3 r = 0.2, the trajectory goes in the non-perturbative 
sector —V" ((fin) — 1 but the IR scaling tends wash out the couplings. 
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4.3. Adding a regulator 

Let us finally discuss how to adapt the results in dimension d ^ 2 when using 
average action. Obviously, one will linearize the flow equation the same way in the 
UV but the loop-integral makes the evaluation of the prefactor of V" more involved 
to establish. The only change is then the definition of /3 r and fi c which are now given 
respectively for Litim, power-law (a p = 1) and exponential (a e = 1) regulators by : 

( d-2\ 1/2 
% = sinc(n—) ft* 

— 1 

= {b-J +1 T(± + 2 ft* (44) 

where the exponential regulator is parametrized by: 

r M = -J£*> (°e> 6 e, Ce) G (M+*) 2 x [1 : oo] (45) 

The behaviour of the /3 C frequency of different regulators with respect to a change 
of dimension is summarized in figure 32. It has been numerically tested that the 




^ 0.8 K 



1.5 2 2.5 
dimension d 



Fig. 32. The /3 C frequency counted in units of the Wegner-Houghton f3™ h frequency (43) against 
the dimension of the spacetime considered for various regulators for the parameters a e = b e = 
c e = 1 and dp = 1, bp = 2, cf ,eq (31), (44) and (45). Note the common crossing at d = 2. 



features obtained by Wegner-Houghton RG can be qualitatively reproduced, namely 
the presence of the instability in the IR for d > 2 and the rebounce driving to the 
gaussian FP in d < 2. 
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5. Discussion 

5.1. Universality in the broken phase for d = 2 

In this paper, we established that at LPA using Wegner-Houghton RG, the sinc- 
Gordon model exhibits a continuous line of non-perturbative IR fixed points in 
d = 2. It was shown that these fixed-points could not be found using Fourier ex- 
pansion, suggesting that the relevance classification around these new fixed-points 
may not be given by the usual simple UV relevant trigonometric functions. It was 
also demonstrated that the fixed points are independent of the bare sine-Gordon 
amplitude and the universality of the /3-functions in the IR has been demonstrated. 

The set of fixed points found in the broken phase allowed us to define the critical 
exponent for the susceptibility in this model 7 ~ 1 for Wegner-Houghton when fitted 
for (i r > 0.95 from figure 23. Obviously, this value of critical exponent is only of 
little help as it was computed near fi r = 1 where the wave-function renormalization 
effect is expected to be large. The next order of the gradient expansion will be 
addressed in another paper 53 . 

5.2. The convexity issue and the Fourier convergence 

As it is well known, the IR behavior of a scalar theory in its broken phase is usually 
driven by a dynamical Maxwell-cut to recover the proper convexity property. Its 
simplest realization is given by the effective potential Vk = —\ ( t> 2 in the concave 
regions. As it has been detailed for both Wegner-Houghton and Average action RG, 
this approach to convexity is usually identified as a tree-level singularity for the 
flow equation. We showed here that in d — 2 the loop-evolution can mimic this 
convexity approach consistently, without sending the flow into a singular region 
where the flow equation does not hold anymore at least in a part of the broken 
phase P r > 0.55, the action being stabilized by IR fixed points. Due to the large 
degenerated sectors of the effective action, the Fourier serie was shown to exhibit 
a very shallow convergence. In d > 2, the Fourier convergence is obviously still a 
problem close to the convexity limit. Using a global resolution not relying on Fourier 
analysis, we demonstrated that the IR flow may run into a genuine singularity. The 
absence of stabilization mechanism can be explained in two ways. 

First, the dimension 2 is a very particular case as it allows to have stable kinks 
in the broken phase. These inhomogeneous configurations contributions to the flow 
may lend some stability to the IR behavior. The lack of proof of stabilization of 
the flow for too small fi r can be seen whether as a too weak numerical accuracy 
whether as the presence of new zero modes : it is known that for low enough f3 r , 
kink bound-states appear in the spectrum. 

Another possible explanation is the non- vanishing scaling dimension in d =/= 2. 
In this case, the dimensionless f3 frequency decreasing with the flow for d > 2, the 
size of the concave regions in the space of the dimensionless field configurations 
is growing as As a result, this enforces the degeneracy of the action in the 
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corresponding regions when the scale is decreased. On the opposite, in d < 2, the 
size of the field configuration space where the action is not convex decreases with the 
flow, and the convexity limit that we could find starting with j3 r « 1 is expected 
to be reppelling. 

5.3. The mass gap 

As widely known from semi-classical considerations, the broken phase of the sine- 
Gordon model exhibits a mass gap due to kinks and kink bound states. The set of 
fixed points found here are clearly incompatible with it as this would mean that 
the dimensionful curvature around the minima is going to zero when integrating all 
the fluctuations. This is the price to pay for the use of LPA : the kink mass which 
is usually recovered by RG to be an invariant of the perturbative flow with wave 
function renormalization crucially depends on the Kosterlitz-Thouless separatrices 
presence. Nevertheless, we saw that even at LPA, one has some traces of the mass 
of the excitations : the convex regions were strongly enhanced in the IR before 
saturating. 

As a result, we can not hope for a global fixed point for the action when includ- 
ing wave-function renormalization non-perturbatively, as the convex regions are 
expected to exhibit a mass gap. 

5.4. The use of average action 

Another aspect of the present study was to clarify the impact of the regulator when 
using average action as its introduction breaks periodicity. It was shown that for 
any regulator, the phase boundary does not change in d = 2. Focusing on the broken 
phase, we exhibited a similar set of fixed points in the IR that have been shown 
to be bare-amplitude independent for various regulators. This feature showed that 
even if the regulator changes quantitatively the broken phase -the critical exponent 
is 7 ~ 0.4 from the fit for (3 r > 0.9 using Litim regulator- it allows to reproduce the 
main features qualitatively. In other dimensions, even if the frequency /3 C is different 
as the one in Wegner-Houghton, one can demonstrate that the results are the same 
provided the use of the relative frequency [3 r . 

In d = 2, as the broken phase appeared to be involved, being the product of a 
fine balance between convex and concave regions of the potential to stabilize the 
potential above the flow singularity, it can also be seen as an ideal testing ground for 
approximate RG schemes, namely proper-time. In 16 > 21 > 40 ) proper-time flows where 
used to study the broken phase of <f> model where it was argued that there was no 
conceptual problem in reaching the limit k = 0. It would be interesting to check if the 
fixed-points are still present, preventing the appearance of a singularity in the flow 
and if the convex region is always well-behaved. Clearly, for the "exponential case" 
of proper-time, that was found to optimize the computation of critical exponents in 
43 , there is little hope to find the fixed points being found here at a finite distance 
of the convexity limit that is in this case rejected at k = 0. 
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As mentionned in the introduction, the average action has been extensively used 
to study Yang-Mills theories deconfinement transition. Due to the presence of the 
regulator breaking gauge invariance, one can use e. g. modified Ward identities to 
wash out its effects 28 . The problem of such process is that the center symmetry 
breakdown that is monitored to study deconfinement is clearly unconstrained by 
such procedure, being a discrete symmetry. Our present finding that the addition 
of a regulator is qualitatively priceless is the case of the sine-Gordon discrete peri- 
odicity symmetry breaking conforts the use of this approach. 

5.5. The phase structure for d ^ 2 

Finally, even simple perturbative considerations allowed us to demonstrate that 
starting with a perturbative sine-Gordon model always drives us in the broken 
phase for d > 2. The absence of a symmetric phase at LPA order is not really good 
news if we are to mimic SU (2) with the sine-Gordon model, as it would correspond 
to the absence of a confining phase. 

This pushes us to conclude that the Haaron gas is more involved than a simple 
perturbative sine-Gordon and that the higher order harmonics dropped for simplic- 
ity in 30 ' 58 have to be taken into account in d > 2 due to the absence of Coleman 
fixed point needed to use universality argument on the first harmonics. In d < 2, 
the symmetry seems to be always preserved, reflecting the fact that d — 2 is the 
lower critical dimension. 

In any case, we saw that the main features in d ^ 2 are driven by the scaling 
dimension of the field so that these conclusions can perfectly be not representative 
of the genuine behavior when wave-function renormalization is taken into account, 
especially in the non-perturbative sectors. 
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Appendix A. Numerical details 

The numerical study of flow equations is in general an involved problem as we 
study a strongly non-linear partial differential equation (PDE) or an integro-partial 
differential equation (IPDE) when the loop-integral is analytically intractable. Many 
studies made the use of the algorithm of 1 using finite differences in a grid of <j>. 
This method has been shown to be very powerful, especially in computing critical 
exponents 11 . 

We made use in the present study of the algorithm developed by Berzins and 
Dew 12 , 13 and 14 known as C° collocation method. The main reason is because 



February 8, 2012 1:37 WSPC/INSTRUCTION FILE sinegordonl 



38 V. Pangon 

it is much faster when dealing with a very fine dependence on <j> as meshes of 10 4 
points can be computed on a laptop within a reasonable amount of time. It has 
been shown to reproduce accurately the analytical results 37 for the Wilson-fixed 
point potential in d = 3 in 54 . The algorithm relies on two different parts. The first 
one, the most important in our case, is the discretization of the field dependence 
using spline representation i.e. Chebyshev piecewise polynomials e. g. 24 and 35 . 
This parameterization exhibits a very fast convergence with respect to two param- 
eters : the degree of the polynomials chosen on each sub-interval and the number 
of sub-intervals into which the field interval is spitted into. 

Due to the finiteness of the range of <f> for which the flow equation is solved, two 
boundary conditions are needed. In the special case of sine-Gordon model, we con- 
strained the first derivative of the curvature to vanish in <p — due to the manifest 
£?2-invariance of the problem. As we solved the flow on an integer number of peri- 
ods =£■, the other boundary has been chosen to follow the same condition due to 
periodicity. Even if periodicity was not explicitly encoded in our numerical solution, 
we found that solving the flow equations over 10 periods were enough to get rid of 
any finite boundary effect. 

The field discretization was dealt on 2000 sub-intervals, each one containing a 
Chebyshev polynomial of degree 3. For low f3 -when the points density is smaller- 
increasing ten times the number of points as been tested to exhibit the same re- 
sults in Wegner-Houghton case and with the power-law regulators. In the case of 
Litim regulator, we tested systematically the numerical stability against an increase 
of number of points and even when the results slightly changed, it never excessed 
10~ 3 in relative error. 

The scale integration is handled by the DASSL routine 18 using backward differ- 
entiation formula that is known to be efficient in estimating the scale derivative in 
stiff cases. 

Finally, the loop-integrals for the power-law regulator are performed using adap- 
tive Gauss-Kronrod rules. Even if the explicit formulation of the flow equation can 
be worked out, cf appendix Appendix C, we chose to perform numerically the in- 
tegration. First because the use of the analytical formulation does not save much 
computation time due to the appearance of non-trivial functions (mostly the arctan) 
but also because we were not able to compute the explicit shape of the flow equation 
for other power-law regulators. As we compared the different flows for b p — 2 and 
bp = 4, we wanted to have the same degree of accuracy in both cases. Let us also 
mention that in the case of the sine-Gordon model, a fit of the loop integral as a 
function of the curvature is not a well suited strategy as this function would have to 
be fitted for the whole range of the possible curvatures. In the concave regions, we 
saw that it gets close to an instability while in the convex regions, a very accurate 
estimate is needed for large curvature values. 
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Appendix B. Callan-Symanzik RG 

Let us briefly present the so-called Callan-Symanzik RG 38 > 7 - 4 . In this scheme, the 
running quantity is not a coarse graining scale but a bare mass and thus the flow 
equation describes how different bare theories are related through a change of the 
mass and reads : 

*fc r M = Tr (pTTP) ) ( R1 ) 

where the functional kernel of the right-hand side have been demonstrated to be 
related to the proper-time flows generators, 39 . In general, this flow equation still 
contains UV (and IR) divergences as no explicit regularization has been made. Let 
us now imagine that we are able to compute the renormalized action for some 
well-chosen value of k. This flow equation is now nothing than the PDE one will 
have to solve to transport the known renormalized result (usually perturbative for 
k large enough) to a region of physical interest (typically k = 0). Obviously, this 
procedure will work provided that a continuous mapping can be made between the 
bare and the renormalized parameters i.e. provided that no phase boundary were 
encountered. 

Let us now apply this equation on the simple ansatz : 



r k [4>] = J d d x [^d^d^ + v k [<f>]] (b.2) 



where z k is field independent. The flow equation for the potential reads : 

that can be derived from the power-law regulator for a p = 1, b p = 1 using average 
action. This integral is straightforward computed in d = 1 as it is well known 
that no divergences occur in Quantum Mechanics. In d = 2, the integral exhibits a 
logarithmic divergence that should be subtracted : 

^=**/»*ot J>»- (ot 



2 z k » \ k 2 

where the term n = was dropped being an infinite field -independent con- 
stant. This equation is nothing but the Wegner- Houghton equation with a field- 
independent wave function renormalization. At this stage, one can remember the 
perturbative proof of 8 to establish the perturbative renormalization group equa- 
tions : starting with a massive sine-Gordon model, all the computations are well 
defined and the vanishing mass limit is taken only at the end. So that this procedure 
can be seen as a perturbative version of Callan-Symanzyk RG. Non-perturbatively, 
we know this procedure does not meet a singularity on the trajectory as the line of 
fixed points we established for Wegner-Houghton holds and protects the evolution. 
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As a result, we believe this provides a non-perturbative proof of the consistency 
of the approach of 8 that has to be checked otherwise only order by order in the 
perturbative expansion. It also shows that no-phase boundary is encountered when 
smoothly switching off the mass term in the massive sine-Gordon model. 

Finally, let us emphasize that the flow equation have been safely regularized 
by the possibility of dropping an infinite contribution to the free-energy that was 
in our ansatz field-independent. If we were to take into account a ficld-dcpendcnt 
wave- function renormalization, it is straightforward to notice that this simple recipe 
does not apply, not even considering the flow equation for Zk- As a result, it may be 
the sign that perturbative renormalization procedure 8 is lost if we were to study 
more involved ansatz. 



Appendix C. Power-law regulator 

While in general, the loop-integral can not be computed in the flow equation of 
average action, it is the case for the power-law regulator for a p = 1 and b p = 2. 
The first point of interest is that it emphasizes the special role of the quantity k 
as the flow equation has two different analytical expressions that are continuously 
matching at V" = k 2 in every dimension. The equation (27) reads in dimension 2 

50. 



kd k V + 2V 

V">2 K 2 



VV" 2 - 4 



Log 



\V"\<2 



\/A - V" 2 



'v" - Vv" 2 - 4' 

V" + VV" 2 -4, 
V" \ 



2arctan 



4 - V" 2 



(C.l) 



Let us emphasize that at V" = 2, at the junction of the two analytical expressions, 
this expression is not only continuous in V" with a common value but also 
infinitely many differentiable. This property of smooth transition between the two 
regions is preserved in the other dimensions as well. In dimension 3, the flow equation 
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reads 



kd k V + 3V- 



V">2 



= -K* 



V2 



<t>V 
V" 



\j v" + Vv^ 1 







\j V" - Vv*' 2 - 4 \J V" + \Fv^ 



\V"\<1 



, -= (C.2) 
VV" + 2 

Finally, the situation in dimension 4 is a little more complicated as the flow equation 
as a superficial divergence in the large momentum limit. However, this divergent 
term is field-independent and can be safely subtracted as performed in appendix 
Appendix B. Unfortunately, after this subtraction, the remaining can not be re- 
casted easily in terms of elementary functions so that the analytical form of the 
flow equation is rather given for V': 

kd k V' + 4V' - 4>V" 

{V" 2 - 4) 3 / 2 

V" - \/V" 2 - 4 N 
V" + \ / V" 2 -4 / 
2K i V" 1 




(4 - V" 2 fl 2 



2arctan 



V" 



a/4 - V" 2 



(C.3) 
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